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By making use of the Langevin dynamics and its generating functional (GF) formula- 
tion the influence of the long-range nature of the interaction on the tendency of the glass 
formation is systematically investigated. In doing so two types of models is considered: 
(i) the non-disordered model with a pure repulsive type of interaction and (ii) the model 
with a randomly distributed strength of interaction (a quenched disordered model). The 
long-ranged potential of interaction is scaled with a number of particles TV in such a way 
as to enable for GF the saddle-point treatment as well as the systematic 1/N - expansion 
around it. We show that the non-disordered model has no glass transition which is in line 
with the mean-field limit of the mode - coupling theory (MCT) predictions. On the other 
hand the model with a long-range interaction which above that has a quenched disorder 
leads to MC - equations which are generic for the p - spin glass model and polymeric 
manifold in a random media. 



I. INTRODUCTION 



The theoretical description of slow dynamics is a crucial point to elucidate the nature of the glass 
transition in structural glass - forming liquids. One of the commonly used approaches, mode-coupling 
theory (MCT), was from the very beginning designed for the supercooled simple liquids Q , i.e. for the 
non - disordered models (as opposed to the models which contain quenched disorder naturally). Later it 
has been proven that MC - equations become exact for a number of spin - glass models [|-|| as well as for 
the polymeric manifold in a random media |^-pd| (i.e. for the models with quenched disorder) provided 
that the number of variables components goes to the infinite. 

Actually applicability of the MC - equations has been substantially extended to the case when the time 
translation invariance and the fluctuation - dissipation theorem do not hold any more ^| . This striking 
similarity between the models with and without quenched disorder suggests that the effective disordered 
potential (e.g., in a supercooled liquid) is in a sense " self - induced" and the difference between such a 
"self - induced disorder" and the quenched disorder might not be crucial |7]j8| . 

In order to provide some insight into self - induced disorder we employed in ref. (IJ] a Feynman 
variational principle (VP) for a set of interacting particles. Indeed it was shown that the VP is capable to 
treat metastable states of the glass - forming system. The main point in ref. |fT2| was that the partition 
function representation in terms of functional integrals is twofold : (a) either as an integral over the local 
density , p(r), or (b) over the conjugated to p(r) field ip(r). It has been assumed that the component 
average free energy , F, (which is only meaningful in the supercooled regime) is equal to the variational 
free energy Fyp. There are at least four strong reasons in favor of this (at first sight not obvious) 
assumption. 



The variational free energy Fyp is a upper bound for the canonical free energy, i.e., F c < F = Fyp 



as it should be since F c = F — TS , where the complexity X > |13|. 



• After implementation of VP the initial problem is reduced to a self - consistent random field 
Ginzburg - Landau model (RFGLM). Then, as was shown previously, the corresponding field ip(r) 
must be upgraded to a replicated field ip a (r), where a = 1, .., n (with the final limit n — > 0) and the 
density field p(r) plays the role of an "external" field. Note that here the density p(r) is Gaussian 
due to the use of the VP. Eventually the correlators of ip and p - fields can be determined self - 
consistently. 

• The resulting replicated partition function for RFGLM has a typical form which may eventually 
lead to the replica symmetry breaking (RSB), structural glass transition and the "self - induced " 
disorder. 
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• Finally, in the case of the long - range interaction the partition function allows the expansion around 
the saddle point , or mean - field (MF), solution. It is possible to show then that the next to the 
mean - field approximation and VP merge and both become exact, i. e. F = F c and the glassy 
phase does not appear. 

Some evidence for this behavior was deduced from the results for the particles on a M- dimensional 
hypersphere Jl4| at large dimensions, M — ► oo. 

The aim of this paper is to face the full dynamical problem for a non - disorded model with a long 
range interaction. Using the expansion around the saddle point solution we derive the full equation of 
motion for the time dependent density - density correlator and show that a "glassy" solution does not 
exist. Conversly, if we add a term describing quenched disorder, by random distribution of the strength 
of the interaction potential, then the resulting equations of motion for two time density correlation and 
response functions fall in the same class as MC - equations which have been widely discussed P-[Tl|] . This 
means that the "self - induced" disorder is not generic for the pure model with the long-range interaction 
and vice versa on addition of a quenched disorder the phase space becomes very rugged resulting in slow 
dynamical processes. 

The paper is organized as follows. In the next section [n] we first introduce the theorotical model 
without quenched disorder . Its dynamics is discussed by using the functional integral technique. The 
saddle point solution yields the mean field dynamcis. Expansions around the saddle point yield one 
loop corrections. The Legendre transformation provides the possibilty of the analysis of the full dynamic 



correlation matrix. In the section [II quenched disorder is introduced by a "random bond model" and 
a Gaussian disorder. The corresponding generating functional (GF) is computed by the self-consistent 
Hartree approximation, which results in a set of coupled Langevin equations, that are solved in their 
asymptotic regimes. More details on the calculations are laid out in the corresponding appendices. 



II. THE MODEL WITHOUT QUENCHED DISORDER 

We start from a simple model system which consists of interacting particles. To do so, let us consider 
a set of N{^> 1) particles in d - dimensional space interacting by a pair potential of the form 

= (£N «p(-«r) 

This is a typical example of a long-range potential with a characteristic length k and a coupling constant 
H/N, The choice of this potential is twofold. It contains a cut off at a k _1 and allows thus to control 
the range of the interaction. Moreover at small scales (r < k _1 ) it consists of a typical power law decay 
with long range character, if < a < 2. Therefore the so chosen potential allows to keep control on 
range and nature of the interaction, which will become essential below. To ensure extensivity of the total 
interaction energy we require that the integral J d d V(r) — O(N ), i.e. it does not depend on the number 



of particles N. As a result we have k cx N a \ The intermolecular potential (2.1) has the form of 
the generalized Kac potential 

V(r) = K d f{nv) , (2.2) 

which has been used for the rigorous treatment of the van-der-Waals theory ]lq ] . In order to provide 
conditions for the expansion around a saddle point, carried out later on (see below), we should require that 
the length must be larger compared to the characteristic size of the system, which scales naturally 
as N x / d ) at N — > oo. As a consequence we find the limits for the range parameter a 

0<a<d. (2.3) 

Below we shall restrict our considerations to the case: d = 3, a = 1, and the strength of the interaction 
fj, > (pure repulsion) without loss of generaliza tion in the main statements we are going to predict. 



Then the Fourier transformation of the potential (2.1) takes especially simple form 



V ^ = (&) • (2 - 4) 

which allows accurate analytic calculations. In the limit N — > oo we have thus k 2 oc TV -1 , but the relevant 
minimum wave vector is k 2 ain oc TV -2 / 3 and thus n 2 can be actually neglected under the integration over 
the whole k - space. As a result we arrive formally at a one - component plasma model (OCP) |lq] where 
the electroneutrality is implicitly provided by a neutralizing background. 
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A. The generating functional method 



In the following we set up the relevant equations of motion for the model system. We restrict ourselves 
to the Langevin dynamics, which can be comfortably formulated in terms of dynamic functionals which 
allows the sys tem atic 1/N - expansion treatment. The Langevin dynamics of N particles interacting via 
the potential (|2.l|) (at d = 3, a = 1 and \i > 0) is described by the equation of motion 



£v«(i 



(p) 



(2.5) 



where niQ and 70 are the mass and the friction coefficient respectively, p — 1,2,..., N and v(r; k) 
cxp(— Kr)/4-Kr. The random force in eq.(|2.5|) is Gaussian with \fi P \t)\ = and the correlator 



/f p) W/j"V)) = 2T l0 5 vm 5 ij 5{t - t') , (2.6) 

where from now on we work in units where the Boltzmann constant ke = 1. 

As was mentioned, it is more convenient to reformulate the Langevin problem ( |2.5| )-( j2~6| ) by using the 
celebrated Martin-Siggia-Rose generating functional (GF) method [0. The method was first applied 
for the (f> 4 - model with the long-range interaction in p8| and for the polymer melt dynamics in |Tg| , |2"c| . 
Despite the fact that the Langevin equation (2.5) is of the second order it is possible to show that the 
Jacobian which appear under transformation to the func tio nal v ariables, still equal to one (see Appendix 
in p3]). After using this technique for the problem (2.5) - (|2.6|), GF takes the form 



N ( N N N 

Z{...}= I f[Dr^(t)Dr^(t)e X plj2A [r^,r^]+ I dt^Y, (t)^ ] v 

J p=l lp=l •* p=lra=l 



(2.7) 



where the action of the free system 



A [rW,fW] = / dt\ T 7o \i&\t)\ + i&\t) 



'dt 2 i 



A Cp)/ 

'dt r ' 



(2.8) 



In the following we are going to transform this functional to collective density variables. By using the 
transformations to the mass density 



rtr) = £*(r-r«(f)) 



(2.9) 



and the longitudinal projection of the response field density 

N 



7T 



r)=E^ P) W V ^( r - r(P) W) 



for the GF one gets 
1 



Z{x«} = J f[Dp a (l)exp^W{ Pa }- ~ j dU2p a {l)U a p{l,2) P0 {2) + J dl Pa (l) Xa (l) 



(2.10) 



(2.11) 



where the summation over the repeated Greek indices is implied. In cq.(2.11) we have introduced 2- 
dimensional field 



Pa(l) 



P(l) 
TT(l) 



(2.12) 



where a = 0, 1 and 1 = (r, t). The "entropy" of the free system is given as usual by 
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W{p, 7r} = log J f[Dr^(t)Cf^(t)exp|^Ao{r (p) ,f (p) } 

N 

p(r,t)-^<5(r-rW(t)) 
7r(r,t) - ir^^VjSir - r W(*)) 



x 5 



x 5 



and [7 a/ g is the 2x2- interaction matrix 



0=1 



EMM) 



Vdri-ral) 
V(|ri-r 2 |) 



(2.13) 



(2.14) 



and Xa(l) is a source field. 

An alternative valuable representation of GF can be obtained through the "functional Fourier trans- 
formation" 



and its inversion 



exp{F{> Q }} = j Dpa(l)exp|w{p a }-» J dlp a {l)^ a {\) 
exp{VF{p Q }}= /" 2ty«(l)expjF{^ a } + i J dip* (l)^a(l) 



(2.15) 



(2.16) 



The substitution of eq.(2.13) into eq.(2.15) leads to the explicit expression for the free-system GF 

N t N 

exp{F{V Q }} = / JjDrW(t)DfW(i)exp^ ^A [rW,fW] 
^ p=i I P =i 

AT „ iV „ 



p—1 p—1 

where ^(1) and (j>(l) are components of the column - variable 

_ ^(i) 



(2.17) 



r=rO)(t) 



(1) 



(2.18) 



By making use (2.16) in (2.11) and after functional integration over p a (l) one gets 



Z{ Xa ,\ a } = f Jj£>Va(l)exp + 

a=Q [ 

dld2[»^ a (l)+Xa(l)] [^ 1 ] Q/3 (l I 2)^/3(2) + X /3(2)]+ f dlVa(l)Aa(l 



(2.19) 



where we have also add a source field A a (l) conjugated to V'a(l)- As a result eqs.( 2.1l| ) and (2.19) 
provide two equival ent r epresentations of GF. For the purpose of expansion around the saddle point we 
use representation ( 2.1 9| ) at A Q (1) = which after the transformation ip a — ► ip a + iXai yields 



Z{x«} = J ]\Dil> a (i)expi-NA fo^JXa] j 



(2.20) 



which is appropriate for a saddle point integration, since the particle number N is large. The action 
hereby is given as 



A& a ;Xa] = \ fdt f d 3 rd 3 rV Q (r,t) [v" 1 ] a0 (r - r'; k)^,*) - ^ log / J] (*)£*<*>(*) 

x exp J £ A [rW , f «] - i £ / dtr« (t) (>> (t)J + z Xct (r^ (t)) ] , (2.21) 
I p=i p=i"' J 



and the interaction matrix 



w Q/3 (r; k) 



1 \ exp(— Kr) 



10/ 47rr 



(2.22) 



Recall that the relation k oc JV is necessary for the validity of the saddle point integration. Moreover 
we have defined the column- vector 



1 



s(p), 



-irj>{t)jdr 



Srf ] {T) 



(2.23) 



for convenience. 



B. The saddle point solution and expansion around the SP 

Minimization of A [ip a ; Xa] with respect to ip a (l) leads to the SP - equations for the mean fields V'a(l) 

(2.24) 

where the average (■■■) SP is calculated by using the cumulant GF 



1> a (T,t) = ~ / d 3 r'v a p{r-r')(pp(r',t)) sp , 



P S p{Va+%} = ^log / JjDrW(t)I3fW(t)ejtpJ £4,[rW 

p=i I p =i 

*E / d<r « } w [^( r(p) w) (r (p) W) 



(2.25) 



The correlation matrix in the random phase approximation (RPA) is defined in such a way 

"<5 (PaCOV 



5 aj g(l, 2) = lim 



N 6x0(2) 



(2.26) 



After linearization of eq. ( p. 25 ) with respect to -0a + *Xq the 2x2 - RPA correlation matrix is easily found 
to coincide with the well known form |l9[ ] 



S a p(l,2) = { F~ l +nv 



(1,2) 



(2.27) 



a,9 



where -C is the interaction matrix (2.22) and F a p is the correlation matrix for the free system F a p(l, 2) 
(Ap a (l)Ap ( g(2)) /JV has the form 



F a/3 (1,2) 



Foo(l,2)F i(l,2) 
F 10 (l,2) 



(2.28) 



In eq.(p8|) Foi(l,2) and Fio(l,2) are response functions whereas Foo(l,2) stands for the correlation 
function. The relation between them is given by the fluctuation dissipation theorem (FDT) which in 
(k, t)- representation has the form 



/3^F 00 (M) =F 01 (k,i)-F 10 (k,t) 



(2.29) 
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It is easy to check that in this case the FDT for the RPA type correlation matrix (2.27) also holds 



- f3-S 00 (k, t) = 5 i(k, t) - 5 10 (k, t) , 



(2.30) 



where f3 — 1/T is the inverse temperature. The corresponding elements of the RPA - matrix ( 2.27 ) are 
of an especially simple form in the Fourier - (k, uj) - representation, namely 



Soo(k,w) = 
5oi(k,w) = 
S 10 (k,uj) = 



F 00 (k,u;) 



[1 + v(k)F w (k, w)] [1 + v(k)F i (k, w) 
F i(k,w) 



l + t>(fc)fbi(k,w) 

F 10 (k,^) 
l + «(fc)Fio(k,«) 



(2.31) 
(2.32) 
(2.33) 



It turns out interesting to recover the wellknown form in the static limit, where we have Sqi (k, lu — > 0) 
/35 s t(k) = (/3F st ) _1 + ju/j -2 and for the correlator S , RPA(k) = S , st (k)/p one gets 



^RPA(k) = 



1 + -p - 



(2.34) 



where we have used F st = p . This expression is completely equivalent to the correlator for the 
OCP-model (see eq.(10.1.7) in EG]) with the direct correlation function c(k) = —fj,(3/k 2 and the De- 

1/2 — 

bye wavenumber = (fl^Po) 



Now let us expand the action (2.21) around SP-solution (2.24) up to the second order with respect to 
the fluctuations tp a (l) — After the functional integration we arrive at the following result for the 

GF 



P{ Xa } =-logZ{ Xa } 

= P S p{^ + Xa}-^Tr[logT Q/3 (l,2)] , 
where T a p(l, 2) is the inverse matrix of the effective interactions p3] 



T a p0;2) = i [^^(l^) + i (Ap Q (l)A P/3 (2)) SP 



(2.35) 



(2.36) 



In eqs. fl2.35 ) - ( |2.36 ) we deliberately keep external field XaO-) nonzero because it to be used in the next 
subsection for the Legendre transformation. 



C. The Legendre transformation 



The functional Legendre transformation is a general way to provide the Dyson equation for the full 
correlation matrix G a p(l,2) p2]. In doing so the irreducible GF, T {(p a (l))} , is defined by the identity 



T{( Pa (l))} + P{( Xa (l))}= / dl<Pa(l))Xa(l) 



By doing functional differentiation of (2.37) one gets 



and 



[G- 1 ] a ,(l,2) = 



S{p«(l)) 

8{ Pa (l))5{pp{2)) 



(2.37) 



(2.38) 



(2.39) 



Taking into account the result in eq.(2.35) we find the following result for GF 



G 



where 



r{(p a (l)» = r SP {( Pa (l))} + ^Tr [logT^(l, 2)] , 



TSP {(Pa(l))} = -P S P {Xa} + / dl ( Pa (l)) • 



(2.40) 



(2.41) 



In eq.( 2.40 ) one should consider X a (l) as a functional of (p Q (l)) given by eq.( 2.38j ). Double differentiation 
of eq.(2.40) leads to an equation of the Dyson form 



[G- 1 ] a M,2)= [S- 1 ] (l,2)-E o/J (l,2), 



(2.42) 



where the RPA - correlation matrix, S a p(l,2), is defined by eqs.( [2.31 ) - ( [2.33 ) and the "self - energy" 
functional S Q| g(l,2) has the form 



E a| a(l, 2) 



-Tr 



27V \S{p a (l))5(pp{2)) 



logT 7(5 (3,4) 



(2.43) 



x==o 



In eq.([2.43|) the "trace" is taken over the variables 3,4 and indices 7, 5. The explicit differentiation in 
eq.( 2.43| ) leads to the result 



E^(l,2) = -— Tr If- 1 



S 2 T 



2N~ 



5{p a (l))6(pp(2)) 



(2.44) 



where T is a short hand notation of the matrix T 7< s(3, 4) and we have took into account that 
<5T Q/3 (1,2)/<5 X7 (3) = (A Pa (l)A Pl3 (2)A Pl (3)) SP /N = at Xa — because the fluctuations are Gaus- 
sian. Further calculation yields 



<5 2 iy(3,4) 
S(p a (l))6( P0 (2)) 



d5d6l(Ap 7 (3)Ap 5 (4)A P( ,(5)Ap x (6)) Sp i? L , (3 (5,2)i? xa (6,l) , 



where 



and the full mean field 



R a p(l,2) = 



ggo(l) 



tf a (l) = -# a (l) + Xa (l) 



(2.45) 



(2.46) 



(2.47) 



The expression for i? a g( l,2) ca n be easily foun d by differentiation of eq.( 2.47 ) with respect to {pp(2)). 
Taking into account eqs. ( g3§ and ( |2~39| ) at x Q ^ one gets 



fl Q/3 (l,2)= [G" 1 ] (l,2)-/i / d4d3^(l,4)S* W7 (4,3)i? 7/3 (3,2) 



or finally 



R af3 (l,2) = / d3< 1 + pvS 



(1,3) 



CK7 



cr 



7/3 



(3,2) 



(2.48) 



(2.49) 



where the h atted variables stands for the corresponding 2x2 matrices. Substitution eqs.( 2.4£ ) and ( 2.45 ) 
in eq.(2.44) yields 



£ Q/3 (1, 2) = - J d3d4AV(3, 4) [G _1 ] ia (3, 1) [G^ 1 ] Sfj (4, 2) 



(2.50) 



where 2x2 vertex - matrix has the form 



K a0 (l, 2) = / d3dM5d6 \ (pv)- 1 + S 



1 + pvS 



(4,3)5 7 t x (3,4,5,6) 



<57 



(5,1W 1 + 



(6,2) 



(2.51) 



X0 
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In eq.(2.51) 5^,^(1,2,3,4) is the 4 - point (response) correlator matrix in RPA 



5^(1,2,3,4) = — (Ap Q (l)Ap^(2)Ap 7 (3)Ap 5 (4)) SP . (2.52) 

ion of 5^(1, 2, 3 ; 
seen as a one- loop diagram (see Fig. 1). 
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The explicit calculation of s (l, 2, 3, 4) is implemented in Appendix B. The vertex - matrix can be 



wvww 

1, a 



WWWV 

2,(3 



FIG. 1. Diagramatic interpretation of the vertex - matrix: the rectangle corresponds to 



the wave 



a^S S^> ^i^'^)' * ne dash line - to the effective interaction matrix [(/i"0) 1 + S 
line - to [i + fi/uS] 

The higher loops contributions which include generally speaking 2m 

S^g 1 ^ 2, . . . , 2m), can be also considered, however the "self - energy" has still the same convolu- 
tion structure : £ = G -1 * K * G -1 . Here the vertex - matrix K a p(l, 2) is calculated in RPA only . That 
is why these contributions basically do not change our results. 



point correlators 



As a result in the (k, u>)- representation the Dyson equation ( 2.42 ) with the "self - energy" functional 
(2.50) reduces to a quadratic one 



G Q7 (k,w) [5- 1 ] 75 (k, w )G^(k,^)-G Q/3 (k, w )+X Q/3 (k, w ) = 0. (2.53) 

The coefficients of the eq. ( [2.53| ) trace the problem back to the free system dynamics which is embodied 
in the correlation matrices F a p(l, 2) and F^) g (l, 2, 3, 4). It is not surprising then that a specification of 



the model for the free system dynamics is necessary, before going to the investigation of eq.(2.53) 



D. Analysis of the equation for the full correlation matrix 

As we have mentioned the explicit solution of eq.( [2.53j ) needs the specification of the free system dy- 
namics. Two simple models are most amenable for the theoretical treatment : free diffusion model 
(FDM) and the relaxation time approximation model (RTAM ) Jll||23| . The latter provide more reason- 
able dynamical information also for short time intervals, At < m^pyQ , where the FDM completely failed 



(e.g., the sum rule does not hold). It turns out that upon calculation of the trace in eq.(2.51) the integral 
is ultraviolet - divergent for FDM and only RTAM leads to the finite result. The matrix elements for 
RTAM has the form 

F 00 (k,u) = ,,f^n D 2 v ( 2 - 54 ) 

F i(k,w ) = — — — — 2.55 

— iuj + k z D — uo z tq 

pF st k 2 D 

Fio(k,w) = ^ — —j- j — i 2 - 56 

ILU + k z D — U! Z To 

where we introduced the diffusion coefficient D = T/70, the characteristic time scale tq = TO0/70, and 
F st = pp for the overall density. At tq — we return to FDM. In the case of RTAM the solution of 



eq.( 2.53 ) for the full correlation matrix reads 



1 + Jl - 4 \-iwT c - uj 2 t q t c + Xst] ^01 (k, uj) 
GoiM = V -^T-. 5 , (2-57) 

2 {-IUJT C - LU Z T Q T C + X st (k)\ 

Gio(k,w)=G i(-k,-w) (2.58) 
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and 



G o(k,w) 



2rl 



l+ V /l-4[- 4 WT c -W 2 T T c + Xst 1 (k)]K l(k,w) 

-iwT c -w 2 r T c +x7 t 1 ( k ) 



^oo(k,w) 



Re 



{ V 1 ~ 4 H wr c - w2r ° Tc + ^/(k)] ifoi(k, w)} 



(2.59) 



The explicit calculation of the matrix K a p(k, ui) (see eg. ( 2.51) is g iven in the Appendix C. The overall 
behavior of the correlation function Goo(k, w) according eq. (2.59) is shown in Fig. 2 (at fj, = 10, (3 = 
0.1, po = 1 and tq = 0.1). It can be seen clearly that there no singularity appears at uj — > 0. 




FIG. 2. The correlation function Goo(k,w) vs rescaled variables u-ro and klg, where 
to = mo/70 and l = (t /^7o) 1/2 

The singularity, however, might be responsible for a glass transition. Instead, the low frequency limit 
of Goo(k, uj) slowly changed with control parameters (which is not shown in Fig. 2 ). That means that 



for the non - disordered model with a general repulsive long ranged potential (2.1) the glass transition 
is not generic. This very important conclusion suggest that for the model with a long range interaction 
the phase space is too smooth to show a glass transition. In order to obtain a glass like transition a 
competing interactions or a quenched disorder should be added. This leads to glassy dynamics, as we 
will show in the next section. 



It is interesting to note, that for the generalized Kac potential eq. (2.2), where /(r) and its Fourier 
transformation are positive definite functions, MCT - memory kernel vanishes at n — > [ p5[ . The 
corresponding argumentation is relegated to the Appendix D. The explanation for this result lies in the 
fact that the "cage effect", which is a cornerstone of MCT, is missing in the MF - limit. 

The "glass transition" which has been studied in rc f. |2q ] for the particles interacted via the Kac 
potential (^^) has a completely different nature. In ref. ]26[] the function /(r) has a step form so that its 
Fourier transform /(k) is negative at some value of k . As a result the system becomes unstable and a 
nonuniform configuration where the particles are grouped into "clumps" shows up. It was found that the 
slow dynamics of the MF - model is associated with these clumps and does not touch a single particle 
motion. Obviously it is different from the conventional glass transition jjjj]. 



III. THE STRUCTURAL MODEL WITH A COMPETING QUENCHED INTERACTIONS 



A. Specification of the model 



In the previous sections we have shown in detail that in the absence of disorder the dynamical spectrum 
changed monotonically with a control parameter and no glassy dynamics can be seen. The natural 
question which arises now is: How will the introduction of competing interactions and / or quenched 
disorder affect the dynamics of the system discussed above? To provide an answer to this question we 
will use already existing models of heteropolymers and their disordered two body interaction |21] Bffl. 
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The use of these models and techniques are here natural, since the behavior of heteropolymers is well 
discussed in the literature. In principle, two practical possibilities exist. 



• The strength of the two- body interaction , fi, in eq.(2.1) is now a random function of all pairs of 
the interacted particles, p pm ("random - bond model"). 

• Each particle carry a single "charge" a p , so that fj, pm = £o + ba p a m and each a p is randomly 
distributed ( "random sequence model" ) . 

It turns out to be sufficient for the purpose of this paper to restrict ourself only with the "random - bond 
model" , where \x pm does not depend from the choice of pairs and has a Gaussian distribution 



P{n pm } oc exp 



(AV» — Co) 2 
2^ 



(3.1) 



The competing long range interactions frustrates the system of particles and the question is whether 
a glass transition exist or not. Normally frustration and frozen disorder is enough for the existence of 
glassy phases. Here the problem is more complicated, since the long range nature of the interaction may 
provide opposite effects. 

The averaging over the quenched disorder in eq.( |2.7| ) (after the substitution fi — > /i pm ) can be carried 
out just in the same way as in ref. [pi 



Similarly, typical two-time dependent terms immediately 
appear. They are also bilinear with respect to the forces of interaction Vju(r) . In order to rationalize 
t hese terms it is convenient to introduce (besides the mass density (2.9) and the response field density 
( 2.10 )) the following collective variables 



JY 



Qo(r, t; r', t') = £ 5{v ~ r*(t))6(r> - r' p (t')) 



P =i 



N 



QiM;r',t') 

Q2{r,t;r',t') 



P =i 

N 



N 



Q 3 (r, t: r>, f ) = " £ ^(t')V i( 5(r' - r* (t'))6(r - r?(t)) 



(3.2) 



0=1 



After the introduction of the 4-dimensional column - fields 



Q a (i;i') 



Qo(i;i') 
Qi(i;i') 
Qa(i;i') 
lQ 3 (1;1') 



(3.3) 



where a = l,2,3,4and 4x4 - matrix 
r ab (l,2,3,4) = 



«(1,3)«(2,4) 
«(1,3)«(2,4) 






the whole expression for GF takes the form 

1 3 
















«(1,4)«(3,2) 




(3.4) 



(Z) av { Xa ,H a } = f f[ f[ Dp a (l)DQ a (l; l') exp | W {p a {l); Q a (l; l')} 
~W j dld2 P<-( 1 )U a p(l,2) Pp (2) + J dl Pa (l) Xa (l) 

dld2d3d4Q a (l;2)r ab (l,2,3,4)Q b (3;4) + J dld2Q a (l;2)H„(l;2) 



AN 2 



(3.5) 
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where the entropy is given by 



W{p a (l);Q a (l;l')} -log f f[ Dr^(t)Dr^ (t) exp {A {r^ , 

J p =i 

Q = 



p a (l)-^rW(lH ri -rW(t)) 



a=l 



JV 



Q a (l; 2) - V pW(l; 2)<5(n - rW(ti))^(r 2 - r«(f a )) 



p=i 



(3.6) 



We had used the column - operators 



ri p) (l) 



iff (<0Vi,i 



pi p) (i;2) 



V 



-irf (f 1 )V i)1 if ! (p) (t 2 )V i , 2 

S (P) 



irr(«a)V A : 



(3.7) 



and the external field, H a (l; 2), co njug ated to Q a (l; 2), has been introduced also. 

The two- point collective fields have a meaning of the dynamical "overlaps" . It is a dynamical 

generalization of the Parisi "overlaps" in a replica space |33|]. For example Qo(l;l') quantify density 
- density a nd Q 2 (l; lQ response - density overlaps respectively between two space-time points. The 
"entropy" fl3.6| ) corresponds to the volume in the dynamic al p hase space when not only fields p Q (l) but 
also overlaps Q a (l; 1') are given. In a sense the "entropy" ( |3.6| ) is again the generalization of the entropy 
for the heteropolymer spanned in a replica space at the given set of "overlaps" p9] . 



B. The saddle point treatment 

Let us introduce the functional F {tp a (l); $ a (l; 1')} by the functional Fourier transformation 

exp{W{p Q (l);Q a (l;l')}} = I JJ II D^ a {l)D<^ a {l; l') exp l')} 

+ i ( dlp«(l)ip a (l)+i f dld2Q.(l;2)$ a (l;2)l . 



(3.8) 



After substitution in eq.(|3.5|) and integration over p a and Q a (l; 2) one gets 



1 3 



(Z) m { Xa ,H a } = f nn^' Q (l)^a(l;l')expjF{^(l);$ a (l;l')} 



N 
N 



dl<Ety a (l) [tT 1 ] (1,2)^(2) 



/" dld2d3d4$ a (l; 2) [r- 1 ] ab (1, 2, 3, 4)$ b (3; 4) 



(3.9) 



where 



F{^ a (l);$ a (l;l')} = log / JjDrW(*)DfW(t)ejq)| ^,4 [rW f«] 

p=i I p=i 

iv . 

-i£ / dir<?>(t) [tf a (rW(t)) (r«(t)) 



JV 



^ / dtdt'pM(t;t') 4> a rW(t);rM(f) 



r( p )(i);r^(i') 



(3.10) 
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In order to ensure the extensivity of the whole effective a ction in eq.( |3.9| ) we put the variance x 2 = x^N 
(so that the variance of the whole streng th fa ctor in eq. scaled as A" 1 / 2 akin to ref. Q) . This 
enable to represent GF in a similar to eq.(2.20) form 



where 



A[ip a ,$ a ;x«,H a } = — I dt I dld2^„(l) [w" 1 ]^ (1,2)^(2) 



\ / dld2d3d4$.(l;2) [r- 1 ] ab (l,2,3,4)$ b (3;4)-llog / f[ Dr^(t)Dr^(t) 

(- iV JV 

xexpi ^^o[r W ,f^]-i^ / dtrW(t) (r (p) (i)) + (r (p) (t) 



p=i 

W r r r 

-*£ / dtdt'p{P\t;t') rW(t);rW(f) 

P =i ^ 



rW(t);r&0 (f') 



The resulting SP - equation reads 



^a(l) = -^ 1 d2v aP (l,2) (p /5 (2)) SP 

1^(1) = -^ 1 d3d4r o6 (l > 2,3,4)(p 6 (2)) SP , 



where the average < • • • >sp is calculated with the GF 



{Z ) av {x a ,H a } = f JJUrW(t)Z»W(t)«pJ £)4>[r« r<*>] 
P =i I p=i 



p=i 



iJ2 / dtdt'pW(t;t') [$ a [r( p »(f);r (p) (t') 



(3.11) 



(3.12) 

(3.13) 
(3.14) 



r (p) (i);r (p) (i')ll [ ■ (3.15) 



Thereby we are left with the GF of a free system which experiences the external mean-fields ip a + i\o 
and W^ + iH a . 



C. The self-consistent Hartree approximation 



In order to calculate GF given by eq.(3.15) we will use the self-consistent Hartree approximation 
(SCHA). For this approximation we replace the real action by an appropriate Gaussian one in such a way 
that all terms which include more then two fields r^\t) or/and f^\t) are written in all possible ways 

as products of pairs of rj it) or f ^ (t) coupled to self-consistent averages of the remaining fields. 

The analogy between SCHA and SP-approximation at N — > oo for the special case when the non- 
quadratic terms in the action are only the functions of the mean-squared displacement d 2 (t — t') = 



Sp=i ( [ r ^(i) — i"( p )(t')] j /N has been proven in ref. ||. In our case the action in eq.(3.15) has a more 
general form. In the Appendix E we show that the SCHA and the next to the saddle point approximation 
(NSPA) merge and both become exact, if the GF with an arbitrary action can be treated by a steepest 
descent approach at N — + oo. 

Let us make the Fourier transformation of the mean-fields 
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1> a (rV>{tj) = I ^3^ Q (k)exp{ikrW(t)} 



a(r (p) (i);rW((')) 



(27T) 

(2tt) 6 



$ a (k\k 2 )exp{ikV p) (t)+ikV p )(i')} 



(3.16) 



(3.17) 



and insert it in eq.( 3.15 ) . Then for eg. ( 3. 15 ) we use the Hartree - type action (see eq.( D24 ). By doing so 
we put for simplicity the expectation value £o = . It is easy to assure oneselves also that the "response 
response overlap" (Qi(l , 1')) = (similar to (ercr) = in rcf. |3l]| ). In the curse of the derivation wc 



have used SP - equation ( |3.14| ) and defined the correlator (or the incoherent scattering function) 



C(k\t;k 2 ,t') = -{Q (k\t;k 2 ,t')) 



as well as the response functions 



G(k\t;k 2 ,t') =--(Q 3 (k 1 ,t;k 2 ,t')) at f < t 
G(k 2 ,t;k 1 ,t') = ^(Q 2 (k 1 ,t;k 2 ,t')) at i' > t , 



(3.18) 

(3.19) 
(3.20) 



where (• • •) stands for the averaging with the Hartree - type of action. After collection of all terms the 
final result (at Xa = and H a = 0) reads then 

r N { N f 



where 



and 



p— l ^ p=i 

dtdt'iY^\t)Y^{t')\(t,t')+ [ dtdt'lY^\t)lY^(t')ri{t,t') \ , 



k 2 \v(k)\ 2 G(k;t,t')C{k;t,t') 



(2tt) 3 

n{t,t') = \xl I ^k 2 \v(k)\ 2 [C(k;t,t')f 



(2tt) s 



(3.21) 



(3.22) 



(3.23) 



In eqs.( |3.2i| ) - ( 3.23 ) we have restricted ourselves to the homogeneous case 

C(k, t; k', t') = (27r) 3 5(k + k')C(k; t, t') 

G(k, t; k', t') = (27r) 3 5(k + k')G(k; t, t') (3.24) 
for the correlation and response function. The equation of motion for the one particle correlator 

W) = j£(rf (*)r?V)) (3-25) 



and the corresponding response function 



m,t') = \Y,{if?{t')rf\t) 

3 = 1 



(3.26) 



(which actually does not depend from the particle index p) can be derived from eq.(B.21) by using the 
standard techniques |H]. The resulting equations read 



d 2 d f 

m °w +l0 di + J rfrA(t ' r) 



V{t,t')- / drX{t,T)V(t,t') 



+ / drr 1 (t,T)g(t',T) = -2T l0 g(t , ,t) 



(3.27) 
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and 



d 2 d 



g(t,t')- I dT\(t,T)g(T,t') 

= -6(t-t') . 



(3.28) 



Eqs.(3.27) - ( |3.28 ) should be supplemented with the initial conditions joG(t + + , t) = —1 and G(t, t) = 0. 
By making use this condition , equipartition (mo/3) Ylj=i ( r j(^)^j(^)) — T, causality G(t,t') = at t < t' 
as well as the condition (1/3) Ylj=i ( r j(t) r j(t)) — one ^ n< ^ from eq.( |3.27 ) the following equation 



1 8 2 
2 m °W 



drX(t,T) 



t t 
dT\(t,T)V(T,t') + [ dTr)(t,r)g(t,T) =2T. (3.29) 



The set of eqs.( 3.27| ) - (3.28) have the sa me s tructure as the Dyson eq.( |2.42 ). After the matrices 
inversions and going to the time domain eqs.( 2.42] ) (in time-translational invariant case) take the form 



m 



m 



dt 2 
dt 2 



7o^ +M(0) 



7o^ +m(0) 



G 01 (t,t') 
Goo{t, t') 



dTX 10 (t-T)G l(T-t')=8(t-t') 

drSio(t - t)G Q0 (t - t') 
J drllnit - t)G w (t - t') = 2T 7o G 10 (t - 



(3.30) 



(3.31) 



where /i(0) = J dt\(t) and RPA - Fourier spectrum 
o 



5*01 



-i^quj — mow 2 + /i(0) 

Eqs.( 3.27 ) - ( 3.28 ) are turned to the Dyson equations ( 3.31 ) - ( ft.3l| ) provided that 

G 0Q (t) = V(t) , G m {t) = -g{t) , 
E 10 (t) = A(t) , En (*)=»?(*). 



We can show 



that the relation 



(3.32) 



(3.33) 



-0— £ii(t) = Sio(t)-E i(t) 
holds, provided that the FDT is satisfied for G a p(t). We then have in addition 

-/^Goo(t) = G i(*)-Gio(*) • 
at 

Bearing eqs.( 3.33| ) in mind the eq.( |3.34 ) takes in our case the form (at t > 0) 

- §- t v(t) = \(t) . 



(3.34) 



(3.35) 



(3.36) 



The validity of the relat io nship ( |3.36 ) can be checked by replacing ( ft. 22 ) and ( |3.23| ) in (3.36). 

The general eqs.( 3.27 ) - ( 3.28 ) are equivalent, mutatis mutandis, to the corresponding equations for the 
p - spin system or a particle in the random potential at the large dimension |7j-|ll|] . The most important 
features of these equations are the glassy dynamical behavior and the universal aging regime. At low 
temperatures the system tries to minimize the energy and each particle (with a number p) tends to 
surround itself with other particles which assure the strength parameter /x pm < 0. On the other hand the 
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long - range interaction tries to support other pairs (ij) corresponding to jiij > 0. As a result the system 
becomes "frustrated" and m any local free energy minima appear. 

In the spirit of ref. |l|,§|||| when t, if — ► oo we have to discriminate between different cases: (i) the 
asymptotic regime when (t — t')/t — ► and (ii) the aging regime when (t — t')/t' — > 0(1). The aging 
regime is much more complicated because the time - translational invariance and FDT are violated. This 
regime has been extensively investigated both theoretically [P, |TT|p4p5| ] and by the computer simulation 
J36| , [37| . In the following we restrict ourselves only to the asymptotic regime, for the sake of clearness and 
simplicity, and since the main features will be already visible. 



D. The asymptotic regime 



This asymptotic regime is characterized by the large time scales, i.e., t, t' — ► oo but keeping the 
difference t = t — t' finite. Under these circumstances we can define 



V &s (t) = lim V(t' + T,t') 

t — -*oo 

£ as (r) = lim g(t' + T,t') 



(3.37) 



Then the equation for the displacement 2? as = 2 [7^(0) — V^t)], response function t/ as (r) and the static 
correlator "P as (0) takes correspondingly the forms 



OT Z OT 



^as(r) - J dr'A as (r - t')V &s {t') - 


00 00 

J air 1 [A as (r + r') - A as (r')] V as {r') - 2 J ' dr' [ V „( T + r') - 7 ?as (r')] &»(t') = 2T (3.38) 



d 2 d , ; 

OT z OT 



r 

Sa S (r) - J dr'A as (r - t')MO = 



(3.39) 



Pas(0) = 



1 



M - A/ a 



00 00 

T - - J drA as X> as (r) - J dT?] as g as (T) 



(3.40) 



where 



M = lim / drX(t, r) 

t^oo J 



(3.41) 



Mas = J dT\ as {r) . 



(3.42) 



However, it is also convenient to define the "anomaly" M = M — Af as O]. The ecjs.(3.38) - (3.40) has 
been analyzed first in the context of polymeric manifold in the random media |9||lC[] and the random - 
phase sine - Gordon model |38| ]. The peculiarity of our model is defined by its memory functions A as (r) 
and r? as (r). 

For example, let as give an explic it exp ression for ?7 as (r). The Gaussian form of the correlator, C(r) = 
exp {— k 2 T> as (T)/2} , leads from eq.( 3.23 )to the result 



f?as(r) 



1 



6 y/V^F) 



(3.43) 
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Usually it is assumed that at the high temperature FDT holds, i.e. 

d 



and 



-13— 2? as (r) =2£ as (r) 



In this case eqs.( 3.38 ) and ( 3.39 ) merge and take a simple form 



mo^T^ + 7o tt + M 
ot 2 - or 



T 

^as(r) - J dr'A as (r - t^V^t') = 2T . 



(3.44) 



(3.45) 



(3.46) 



It turns out |||T(^[37j that the solution which satisfies the FD T is on ly stable above a critical temperature 
T c . For the stability analysis it is convenient to represent eq.(3.46) in the form 



m 7T^ + 7o 7T +M + M as (r) 



V as (r) ~ I At' [77 as (r - r') - ^(r)] ^pDjf) = 2T , 



(3.47) 



where 



oo 

M as (r) = J rfr'A as (r') 



For t — > oo the stability condition which comes out of eq.(3.47) reads 

[M + M as (r)] P as (r) < 2T . 
Then the stationary value of the displacement T> as (r — > oo) = go reads 

2T 



Qo 



AI 



By taking into account eqs.( p.44| ) and ( p.45[) the stability condition becomes 

D(q,T)>0 

for < q < qo, where 



D(q,T) 



\t) 12^ q a 



Xo 
T ) 12 



(3.48) 



(3.49) 



(3.50) 



(3.51) 



(3.52) 



The critical values <7 C and T c at which the condition (3.51) first becomes violated is defined by equations 

(3.53) 



D(q c ,T c )=0 
D'(q c ,T c ) = 

Consequently, eqs. ( 3.53| ) have the simple solution 

a/^90 



24 



and q c = q Q 



(3.54) 



The Fig. 3 shows the behavior of D(q,T) in the vicinity of the critical point. It can be seen that the 
minimum, q m < qo, at which D(q,T) < appears continuously, i.e. the instability of the FDT solution 
shows up as a 2 nd order phase transition. This is analogous to the dynamics of polymeric manifolds in a 
medium with the long range correlation in disorder ]To| . In particular, if "anomaly" M — > then qg — > oo 
and T c — > oo, so in this case the FDT solution is unstable for any finite temperature. 



16 



D(q,T) 
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1000 1050 

q 



FIG. 3. D(q,T) vs g at xo = 0.1, go = 10 3 for different temperatures: (i) full line corre- 
sponds to T = T c = 0.1528; (ii) dashed line T = 0.1535 ; (iii) dot - dashed line T = 0.151 

Let us consider the dynamics at the temperature slightly above the critical point: T = T c (l + e), where 
0<e<l. For large r the decomposition 



2\ s (r) = q - f(r) , 



(3.55) 



is possible, where f(r) -C go- The substitution of this decomposition into eq.(3.47) and the expansion up 
to the second order with respect to /(r) yields 



r 

eqof(r) + ~ [/(r)] 2 + \ J dr' [f(r r') f(r)] A/( T ') = 



(3.56) 



Following ref. |l| let us make the Laplace transformation C[f(r)} = f(z) and introduce the scaling 
functions, 4>{z) or 4>{t) , in a such way 



/(*) 



(z) or /(r) = c s 4>(t) 



(3.57) 



where i — z/uj e and f = tuj e . If c e — e and ui e = loqS 1 / ' then one can write eq.(3.56) in the form 

3 



qo$(z)- -C{$ 2 (r)}(z) + z$ 2 (z) = 



(3.58) 



(see eq.(2.68b) from ref. §). 

In the critical regi me z ^> 1 (or f < 1) the solution of eq.( 3.5S ) has a form <j)(f) cx f _a . In this limit 
the first term in eq.( |3.58| ) is dropped out and the exponent is defined by the equation 



r 2 (l-a) _ 3 
T(l - 2a) ~ 4 



(3.59) 



The solution of eq. ( |3.59 ) gives a = 0.30465. In the opposite limit z <C 1 (or f > 1) the last term 
in eq.( |3.58 ) can be neglected. In this case the solution has the form 4>{t) cx A e r~ a exp{— Af }, where 
A e = 8eg r(l - a)2( 1 " 2a )/3r(l - 2a) A 2 . As a result the overall scaling reads 



2\ s (r) = 



9o 



A 



-, at Lo e T <c 1 



i6 ra exp{— A(w e r)}, at oj £ t 3> 1 



(3.60) 



where A is some constant. 

At T < T c the FDT is violated for the large time separation r and the aging regime is arising. It should 
be mentioned that the asymptotic regime cannot be decoupled from the aging one fHHJ. In actual fact, 
the "anomaly" M in the asymptotic eqs. (3.35) - ( 3.40| ) strictly speaking can be calculated only from 
the aging regime. Because of the distinct aim in this paper we are not going to discuss the aging regime 
here expecting to return to it in a later publication. 
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IV. CONCLUSION 



In the present paper we have considered the dynamics of two models with the long range repulsive 
interaction. The interaction potential was designed in a way to enable the saddle point treatment, as 
well as a fluctuation expansion. 



For the pure model we have derived eq.(2.53) for the full correlation matrix G Q| a(k, u>) in the one-loop 
approximation, which has an the explicit solution (see eqs.(2.57) - ( p. 59 ) ). This solution has a "boring 
behavior" at u> — ► which manifests the absence of the glass dynamics. The physical background of this 
stems from the fact that the potential is much too soft and the "cage effect" is completely missing. 

This conclusion is in accordance with the interacting particles statistical thermodynamics analysis, 
which was given in ref. Jl3|. It was shown there that for the infinite range interaction potential, which 
allows a well defined saddle point treatment, the glassy phase is simply suppressed. 

On the other hand, the same model but with a randomly distributed strength of interaction (the 
"random - bond model" ) leads to the continuous glass transition. This type of transition is also the 
case for the polymeric manifolds in the disordered medium with long range correlation [ [Tof as well as for 
the p - spin interaction spin-glass model at the large external field |3||(|. It would be also interesting to 
investigate the more realistic "random sequence model" when each particle carry a random "charge" . 

Qualitatively the same glassy behavior have been found in the pure spin models with the deterministic 
but very rapidly oscillating coupling between variables |39|-f42]] . It was assumed that the effective quenched 
disorder is "self - induced" [p 39 1 . This means that because of the slow dynamics some degrees of freedom 
freeze and play the role of the effectively quenched disorder. 

As a conclusion, the glass transition in the pure systems of the interacting particles, where the disorder 
is actually "self - induced" , goes beyond the mean- field level . This appears to hard to implement in 
the present context, because it implies the consideration of the short range interaction potential as well 
as activated processes. 
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APPENDIX A: CALCULATION OF THE 4-POINT RPA - CORRELATION MATRIX 



In the full analogy with eq.(2.26) the expression for 5^ 7<5 (1, 2, 3, 4) reads 



lim 



o [NH Xfi (2)5 Xl (3)5 X 6 (4) 



SP 



(Al) 



The expansion of the_Jo a (l)) SP up to the 3 d order with respect to the mean field ip a + %Xa can be easily 
obtained from eq. (|2.25 ) 



(Pa(l)> S p = (pM))o + j d2(A Pa (l)A Pfj (2)) [xp(2) - i^(2)] 

d2d3 (Ap Q (l)A P/3 (2)Ap 7 (3)) [ X p(2) - %(2)] [* 7 (3) - ^(3)] 
d2d3d4(Ap a (l)Ap /J (2)Ap 7 (3)Ap,(4)) 



x [ X p(2)-ifa(2)] [x 7 (3)-# 7 (3)] [ Xs (4)-iM4)] 



(A2) 



By using the SP-equation ( 2.24 ) and after threefold differentiation with respect to Xc*(l) ( see eq.( Al)) 
we find 



5^(1,2,3,4) = FW (1,2,3,4) 



1 + nvF 



(2) 



(1,1) 



1 + pvF 



(2) 



(2,2) 



00 
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(2) 



1 + fj,vF 



where the 4-point free system correlation matrix 



(3,3) 



1 + nvF 



(2) 



(4,4) 



F, 



(4) 



1 



^(1,2,3,4) = (Ap a (l)A P/3 (2)Ap 7 (3)Ap A -(4)) 



(A3) 



(A4) 



In eq.(|A3j) we imply the sum mati on (integration) over the barred indices (barred space - time 
variables). When deriving eq.(A3) we have also kept in mind that (Ap Q ,(l)Ap / g(2)) oc N and 

(Ap Q (l)Ap /3 (2)A / 9 7 (3)Ap5(4)) oc N 2 , etc. The fact that the matrix F^^l, 2, 3, 4) is symmetrical 
with respect to simultaneous permutations of Greek indices and space-time arguments as well as eq.(2.27) 
have been used. 



It is easy to show that F^g s (l, 2, 3, 4) is factorized 



flL (1,2,3,4) = ^(1,2)^(3,4) +Fi%l )3 )^(2,3) + ^(1,4)F^(2,3) 



,(2) 



(2), 



p(2), 



a/3 



7<5 



Pi 



(A5) 



On the other side it is instructive to check that even in this case S^} g (l, 2,3,4) can not be factorized 



APPENDIX B: CALCULATION OF THE VERTEX-MATRIX K a p{l,2) 



The substitution of the eg. (A3) into eq.( 2.51 ) after the straightforward algebra yields 



K af) {l,2) = \ (2//) 



-1 A-l 



F 



/i -1 *) -1 +F 



(2,1) 



{^(l,2)F 7 j(3,4)+^ a7 (l,3)^2,4)+^ 



{[l + ^vF] l \ (3,1) 



1 + 2[ivF 



(4,2) , 



(Bl) 



5/3 



where as before for the repeated barred indices (variables) the summation (integration) is implied. For 
the time - space - translational invariant case the respective Fourier transformation leads to the result: 



K, 



a p(k,Lu) = |/F 7 j(k,w) 



(2ft) 1 t>- 1 +F 



(2^*) 1 v~ 1 +F 



-1 


/T 1 *) -1 +F 




• (-k,-w)F a? (-k,-w)%(k,w) 








/3a 


-1 


H^v- 1 +F 




• (k,w)F aJ (k,w)F^(-k,-a;)} 








/3a J 



{[l + 2/ifiFl H (-k, -w) { [1 + 2/ivF] H (k,w) , 



where the trace 



d 3 qdu 

~(2^y 



1 + /i^F^lT 1 



(q,^) 



(B2) 



(B3) 



With the correlation matrix F given by eqs.( [2.5q ), by doing integration over u> one can check that the 
trace 1 = 0. This gives finally 



(B4) 



where 



L a/3 (k,w) 



1 + 2[ivF 



(k,w)F^(k,w) 



«7 

l + {2^)- 1 F- 1 v' 1 
l 



1 + /i" 1 F" 1 ti" 1 



(q,^) 



/35 



1 + 2/xt)F 



(k,w) 



(B5) 
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APPENDIX C: THE MCT FOR THE GENERALIZED KAC POTENTIAL 



In this case the direct correlation function c(r) = — 0V(r) and its Fourier transformation takes the 
scaling form 

c(k) = -/?/ (^j . (CI) 

Let us insert this expression to the MCT - memory kernel (see eq.(3.32) in JO). It is reasonable then 
to rescale the integration variables in the memory kernel , k — > ftk, p — > up as well as to put for the 
external wave vector q = «qo , where qo is some reference wave vector. The last scaling means that in 
the MF - limit an experiment probes a very small wave vector : q — > 0. The resulting scaling of the 
memory kernel , m(q, i), reads 

m( Kqo ,t) = K d S(q ) ^ |^^(k + p-q ) 

{ e L (q)/3[k/(k)+p/(p)]} 2 ~ ~ 

x -i '— 5(k, t) S(p, t) , (C2) 

% 

where we have took into account the scaling form of the correlator: 5(k, t) = S(\t/n; t). Thus we finally 
arrive at the scaling m(Kq ,t) cx n d — > and the glass transition dies out. 



APPENDIX D: THE ANALOGY BETWEEN SCHA AND NSPA 



Let us prove that SCHA becomes exact for GF given by (3.15) in the limit N — > oo. We will consider 
even a more general GF 

Z{ Xa }= /il II Dx « ] W ex p|-^E / dU2x { i\l)A otP {l,2)xf{2) 

J p=la=0,l I p=l ^ 

+ E M/ K P) ]+E J dlx^(l) Xa (l) \, (Dl) 



p=1 p=i 
where we have used the shorthand notations 



■' :; " (M l^p 



and "1" embraces Cartesian indices as well as time variable: 1 = {i,j,k;t}. In cq.(Dl) 
arbitrary functional of . 



» 



(D2) 



is an 



Instead of the exact action functional in eq.(Dl) we consider now the trial one which has a Gaussian 
form 



A' 



1 



P =i 



L P) (1) =E o / ^^(1)^(1,2)^(2) 



dld2 ^(1)^(1,2)4^(2)- / dlL a (l)a#>(l) 



(D3) 



Let us look for the "best" coefficients r Q/ g(l, 2) and L a (l) in a sense that the exact "free energy" F [x a ] 
— log Z {xa} tends to the trial one F [ Xa ] = — log J ]1 Dx^ exp{— S[a;^]}, i- e. 



F \Xa] — » ^0 [Xa] 



(D4) 



and both becomes exact at TV — > oo. 

We can show that the property ( |D3| ) is satisfied by T a p and L Q which are obtained by extremization 
of the functional 
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N 



$>{T al 3,L a } =-log / J] II Dv&O-) e*p\-S 

J p=la=0,l I 



N ( . 

+ J2\ dld2r^(l,2) (^(1)^(2) 
p=i I J 

+ jdl [L a (l)-x a ](x%\l))-{w[x%){l) 



(D5) 



where we use the notations 



/n; =i n a =o,iM p) (l) ...exp{-5 


'JpY 


} 


/n^ 1 n a =o,i^ ) ( 1 ) ex p{-' s 


'(p)~ 




> 



(D6) 



The extremization conditions reads 



<sr Q7 (i,2) 



SL a (l) 



$ = 



$ = 



(D7) 



The variations in eqs.( p6| ) can be done directly. During the calculation th e ge neralized Wick's theorem 
pg| should be also taken into account. Namely, because the averaging (D5) is simply the Gaussian 
integral the Wick's theorem yields 



+ / d2 (A^)(1)A4 P) (2) 



W 



M 



(D8) 



where Ax^(l) = Xa'(l) — (x { a'(l) ) . After the straightforward calculation we find 



Mi 



Mi 



r a7 (i,2) = - 



2 \8x [ S\l)6xf{2) 



W 



riP) 



and 



L a (l) = j d2[i4 a/9 (l,2)-r a/3 (l,2)] (a#>(2)) 
Then equations for the two moments take the form 

' 5 



(*«(!)) =/d2 [A" 1 ] ^(1,2) 



and 



A^)(1)A4 P) (2)\ = 



,foW(2) 



A-2T 



W 



,(p) 



X/»(2) 



(1,2), 



(D9) 



(D10) 



(Dll) 



(D12) 



a/3 



where A and f stands for the corresponding 2x2- matrices. 



On the other side the saddle point (SP) treatment of eq.(Dl) at N — > oo yields 



d2 A a/3 (l,2) xf{2) + 



SW 
Sxi p \l) 



+ Xa(l) = 



(D13) 
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and 



A*(f)(l)A4 p) (2) 



sp 



A — 2B 



(1,2), 



a/3 



where 



B a p(l,2) 



8 2 W 



2 Sx { l\l)dxf\2) 



(D14) 



(D15) 



and ir^(l) stands for the field in SP. 

In order to show the analogy between eqs.(Dll)-( D12| ) and eqs.(D13)-(D14) let us make the functional 
Fourier transformation 



exp{x[ 2/ (f)(l)]} = J Dx^(l)explw[x^(l)\ -i j dlx^(l)y^(l] 
and its inversion 



exp ■ 



■{w[x^(l)]}= j DyW(l)exp|if [^(1)] +i j dl x^(l)y^(l) 
Then eqs.( D13 ) - ( D14 ) can be written as 

*W(1)= / d2[A- 1 ] a0 (l,2) \ i (y^(2)) sp+X ,(2) 



and 



As« (1)^(2) 



SP 



i+|| (AyAy) 



SP 



(1,2) 



a/3 



where the correlation matrix 



(AyAy) SP \\^(A y ^(l)Ay^(2) 



SP 



and 



(■■■)sp = 



J ... exp [y<f>] + i / dl 5><f > (l)yi p) (1)} 



y (p) 



At — + oo by making use eq.( |D17| ) one can immediately see that 

/ S 



fdl x& (1)^(1)} 



W 



*(«l)) sp 



and 



rf(l)54 p) (2) 



Ay^\l)Ay^{2) 



sp 



(D16) 



(D17) 



(D18) 



(D19) 



(D20) 



(D21) 



(D22) 



(D23) 



and SCHA exactly corresponds to NSPA . For the case which was treated in Sec.III C (x^\l)) = 
and the Hartree - type action ( |D3| ) cast the form 

N 



S 



4 P) (1)| =EU I dld2x^(l)A a0 (l,2)xf(2) 



0=1 



2 

dld2 



5x ( S\l)5xf{2) 



W 



x ( S\l)xf{2) 



(D24) 
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FIG. 4. Diagramatic interpretation of the vertex - matrix: the rectangle corresponds to 5^.^(1,2,3,4); the 
dash line - to the effective interaction matrix [(pit)) -1 + S 1 ] ; the wave line - to [l + fivS^ 
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FIG. 6. D(q,T) vs q at xo = 0.1, qo = 10 3 for different temperatures: (i) full line corresponds 
T = T c — 0.1528; (ii) dashed line T = 0.1535 ; (iii) dot - dashed line T = 0.151 
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